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ABSTRACT 

This  paper  describes  a  theoretical  investigation  of  the  effect  of  mass  transfer 
into  the  hypersonic  low  Reynolds  number  viscous  layer  upon  the  heat  transfer  rate 
and  skin  friction,  and  represents  an  extension  of  an  earlier  study  by  the  present 
authors.  New  results  have  been  obtained  and  correlated  as  a  function  of  the  Reynolds 
number.  The  range  of  applicability  of  the  present  fluid  dynamic  model  is  also  dis¬ 
cussed  and  comparisons  are  drawn  with  earlier  work,  both  theoretical  and  experimental. 
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mass  fraction  of  species  i 
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chemical  source  term,  net  mass  rate  of  production  of  species  i  per 
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Subscripts 


density  ratio  across  shod:  wave 

transformed  coordinate 
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wall,  surface  of  vehicle 


upstream,  ambient  conditions 
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Le 
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skin  friction  coefficient 


binary  Lewis  number 
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I.  INTRODUCTION 


For  the  controlled  flight  of  slender  ballistic  re-entry  vehicles  and  hypersonic 
lifting  re-entry  vehicles  at  high  altitudes,  it  is  mandatory  that  adequate  information 
be  available  for  the  prediction  and  control  of  the  exchange  of  mass,  momentum  and 
energy  between  the  surface  of  the  vehicle  and  its  environment.  Hypersonic  flight  in 
the  Earth's  atmosphere  from  the  outermost  regions  down  to  the  surface  can,  in 
principle,  be  treated  and  solved  by  utilizing  kinetic  theory  and/or  continuum  me¬ 
chanics.  However,  the  problem  is  greatly  simplified,  if  we  suitably  divide  the 
entire  regime  into  sub-regimes  in  which  different  aerothermodynamic  phenomena 
dominate. 

In  Figure  1,  the  high  altitude  hypersonic  flight  regime  is  shown  divided  into 
sub-regimes  most  suitable  to  the  present  discussion  (Ref.  1).  Typical  ballistic  and 
lifting  re-entry  trajectories  of  hypersonic  vehicles  are  superimposed.  Note  that 
Figure  1  is  typical  for  a  nominal  nose  radius  equal  to  one  tenth  of  a  foot;  for  smaller 
nose  radii,  all  of  the  low  density  effects  are  shifted  to  lower  altitudes.  Descriptive 
stagnation  region  velocity  profiles  are  shown  to  the  right  of  the  respective  flight 
regimes.  These  regimes  are  discussed  below. 

As  is  well  known,  hypersonic  flight  at  low  altitudes  is  associated  with  high 
Reynolds  numbers  and  a  flow  field  defined  by  a  very  thin  shock  wave  followed  by 
an  influscid  region  separated  from  the  body  by  a  thin  boundary  layer.  This  is  the 
ordinary  hypersonic  continuum  regime,  and  the  whole  flow  field  is  known  as  the  shock 
layer.  Conditions  directly  behind  the  shock  wave  are  determined  using  the  Rankine- 
Hugoniot  relations  and  the  influscid  flow  is  determined  by  using  one  of  the  available 
hypersonic  shock  layer  approximations  (see,  for  example  Ref.  2);  e.g. ,  modified 
Newtonian,  or  stream  tube,  in  order  to  obtain  outer  boundary  conditions  for  the 
boundary  layer. 

For  suborbital  flight  speeds,  the  radiative  heat  transfer  from  the  shock  layer 
is  very  small  compared  with  the  aerodynamic  heat  transfer  and  hence  the  dominant 
energy,  momentum  and  mass  transfer  processes  can  all  be  treated  by  studying  the 
phenomena  within  the  boundary  layer  adjacent  to  the  surface.  This  procedure  has 
yielded  many  useful  solutions  to  the  thermal  protection  problem,  (e.g. ,  Refs.  3-6). 
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At  a  somewhat  higher  altitude  (lower  Reynolds  number)  the  fluscid  effects  extend 
throughout  a  major  portion  of  the  shock  layer  and  the  flow  may  no  longer  easily  be 
separated  into  three  distinct  regions.  Convention  has  now  termed  this  flight  region  the 
hypersonic  low  Reynolds  number  flow  regime  which  may  be  further  divided  as  follows: 

1)  Viscous  layer  subregime  -  It  is  assumed  that  the  classical  conserva¬ 
tion  laws  of  macroscopic  flow  arc  applicable  and  the  shock  remains  sufficiently  thin  to 
be  treated  as  a  discontinuity  followed  by  at  least  a  small  influscid  region. 

2)  Merged  viscous  flow  subregime  -  At  sufficiently  low  Reynolds  numbers 
the  thickened  shock  merges  with  the  fully  fluscid  flow  behind  it,  and  the  entire  flow 
field  from  the  free  stream  to  the  body  surface  should  be  treated  as  a  single  problem. 

At  extreme  altitudes,  (the  free  molecule  flow  regime)  it  is  usually  permissible 
to  neglect  collisions  between  incident  and  reflected  molecules  in  studying  interactions 
between  the  vehicle  and  its  environment.  However,  it  is  vitally  important  that  the  state 
of  the  undisturbed  gas  (Ref.  7),  and  the  surface-particle  interactions  (Refs.  8-10),  be 
fully  taken  into  account.  With  decreasing  altitude,  particles  reflected  from  the  surface 
begin  interacting  with  the  incident  particles.  It  becomes  necessary  to  utilize  kinetic 
theory,  e.g. ,  the  Boltzmann  equation  (Refs.  11-14)  in  order  to  obtain  solutions  in  the 
near  free  molecule  flow  regime. 

Between  the  low  Reynolds  number  and  near  free  molecule  flow  regimes,  a  transi¬ 
tion  regime  exists  in  which  the  shock  wave  is  formed  requiring  a  more  rigorous  kinetic 
theory  approach,  (Ref.  15)  possibly  coupled  with  a  macroscopic  analysis  in  which  higher 
order  terms  are  introduced.  It  should  be  noted  that  the  uncertainty  in  the  extent  of  this 
region  is  indicated  by  the  lightly  shaded  area  plus  the  broken  line  in  the  velocity  profile. 

It  is  seen  from  Figure  1  that  knowledge  about  the  hypersonic  low  Reynolds  number 
flow  regime  is  important  both  for  slender  ballistic  entry  vehicles  and  lifting  re-entry 
vehicles.  The  earliest  studies,  restricted  to  order  of  magnitude  analyses,  are  sum¬ 
marized  by  Hayes  and  Probstein  (Ref.  2)  where  the  entire  hypersonic  flight  regime  is 
divided  into  seven  subregimes  from  continuum  to  free  molecule  flow,  including  the 
boundary  layer  regime,  vorticity  interaction  regime,  viscous  layer  regime,  incipient 
merged  layer  regime,  fully  merged  layer  and  transitional  layer  regimes,  first  order 
collision  theory  regime  and  the  free  molecule  flow  regime.  This  classification  scheme 
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has  been  the  basis  of  most  subsequent  studies. 

In  Hoshizaki's  study  (Ref.  16)  of  the  effect  of  bluntness -induced  vorticity  on 
the  boundary  layer,  it  was  concluded  that  the  inviscid  momentum  equations  evaluated 
behind  the  shock  are  equivalent  to  using  the  vorticity  as  a  boundary  condition,  and 
thus,  the  complete  momentum  equations  incorporate  the  vorticity  boundary  conditions. 
Increases  of  the  order  of  ten  percent  in  heat  transfer  and  skin  friction  were  obtained 
when  the  vorticity  terms  were  retained  in  the  Navier-Stokes  equations,  as  the  Reynolds 
number  was  decreased  to  a  minimum  valid  Reynolds  number.  However,  these  terms 
are  retained  in  an  appropriate  low  Reynolds  number  order  of  magnitude  analysis. 

In  Refs.  17-20,  similar  results  were  obtained  by  appealing  to  the  low  Reynolds  number 
nature  of  the  flow  without  special  recourse  to  vorticity  except  as  a  modification  of  the 
pressure  term. 

In  Refs.  1G-19,  the  density  was  assumed  constant  through  the  shock  layer. 

Since  it  can  easily  be  shown  that  the  variation  of  density  through  this  shock  layer  is 
of  the  same  order  of  magnitude  as  that  across  the  shock  wave,  it  would  appear  that  the 
assumption  of  constant  density  is  not  realistic.  In  a  more  recent  study  (Ref.  20),  the 
density  variation  was  taken  into  account,  and  it  was  shown  that  as  the  Reynolds  number 
decreases,  the  ratio  of  heat  transfer  calculated  with  density  variation  to  the  heat 
transfer  calculated  with  constant  density  increases.  Clearly,  then,  the  variation  of 
density  across  the  viscous  layer  should  be  included  in  a  theoretical  model. 

In  all  of  the  above  studies  however,  solutions  were  obtained  for  the  flow  of  a 
viscous  non-dissociating  perfect  gas.  Hence,  these  studies  ignored  the  important 
physicochemical  phenomenon  of  dissociation,  which  decreases  the  molecular  weight 
of  the  gas  and  alters  the  relationship  between  density,  temperature  and  pressure. 

Also  not  considered  was  the  alteration  in  the  flow  produced  by  mass  transfer  from  the 
surface,  which  was  subsequently  shown  to  have  a  considerable  effect  upon  the  structure 
of  the  viscous  layer  and  other  dependent  variables,  (Ref.  1). 

In  Ref.  1.7,  Probstein  and  Kemp  extended  their  analysis  to  the  merged  viscous 
layer  problem  by  matching  a  simplified  perfect  gas  shock  structure  solution  with  their 
constant  density  perfect  gas  viscous  layer  solution.  More  recently,  Levinsky  and 
Yoshihara,  Ref.  21,  also  studied  the  merged  viscous  layer  but  included  variable  density. 
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Qualitatively,  the  latter  results  appear  credible.  However,  in  addition  to  the 
remarks  made  previously  about  the  non-dissociating  perfect  gas  model  in  the  vis¬ 
cous  layer  solution,  its  use  in  the  study  of  hypersonic  shock  wave  structure  in  air 
is  even  more  unrealistic,  since  not  only  is  the  gas  within  a  hypersonic  shock  wave 
not  in  chemical  equilibrium,  it  is  not  necessarily  in  thermodynamic  equilibrium 
(Ref.  22)  i.  e. ,  the  relaxation  times  for  the  translation  degrees  of  freedom  are 
much  shorter  than  the  other  degrees  of  freedom,  (e.g.,  rotational,  vibrational). 

In  particular,  the  available  kinetic  energy  of  the  molecules  is  redistributed  among 
many  different  internal  modes  of  energy  storage  which  absorb  energy  at  different 
rates,  but  nevertheless  in  parallel  processes  and  not  in  scries.  It  is  stressed  that 
the  dissociation  and  ionization  processes  which  absorb  energy  and  act  to  reduce  the 
gas  temperature,  will  begin  before  thermal  equilibration  is  complete. 

Cheng  (Ref.  23)  examined  the  merged  viscous  layer  without  recourse  to  the 
weak  approximation  of  a  perfect  gas  shock  wave  structure.  Instead  as  suggested  by 
Sedov,  et.  al.  (Ref.  24)  the  Rankine-Hugoniot  shock  relations  were  modified  by  assum 
ing  first  order  corrections  to  account  for  heat  conduction  and  the  normal  component 
of  momentum  transfer  due  to  the  viscous  effects  in  the  shock  wave.  With  this  ap¬ 
proach,  he  has  apparently  been  able  to  extend  the  perfect  gas  low  Reynolds  number 
viscous  layer  analysis  out  to  almost  the  perfect  gas  near  free  molecule  results  of 
Willis  (Ref.  14). 

The  present  analysis  has  been  carried  out  in  the  regime  defined  in  Figure  1 

as  the  hypersonic  low  Reynolds  number  viscous  layer  regime.  In  this  regime,  the 

viscous  layer  thickness  6  is  significant  compared  to  the  shock  detachment  distance 

6  ,  but  6/6  <1,  and  the  shock  wave  thickness  is  assumed  to  be  sufficiently  small 

s  s 

so  that  it  can  be  assumed  that  the  flow  is  in  thermochemical  equilibrium  behind  the 
shock,  before  entering  the  viscous  layer. 

The  governing  equations  utilized  in  the  analysis  are  the  equations  of  change 
for  the  flow  of  a  compressible  chemically  reacting  gas.  Included  are  the  conserva¬ 
tion  of  mass,  momentum  and  energy.  In  order  to  simplify  the  system  of  equations, 
the  diffusion  equations  were  uncoupled  by  assuming  a  Lewis  number  of  un*ty.  The 
boundary  conditions  are  split;  some  are  specified  at  the  surface  and  the  others  behind 
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the  shock  wave. 

To  evaluate  the  remaining  transport  coefficients,  the  Prandtl  number  was 
assumed  to  be  constant  at  a  value  of  0.71  and  Sutherland's  law  was  utilized  for  calcu¬ 
lating  the  viscosity  of  dissociated  air. 

Even  after  expanding  the  simplified  version  of  the  governing  equations  into 
the  usual  body-oriented  coordinates,  it  is  necessary  to  make  further  simplifying 
assumptions  in  order  to  reduce  the  system  of  equations  to  a  more  tractable  form. 

It  is  further  assumed  that  the  shock  layer  thickness  ft  is  small  compared  with  a 
characteristic  dimension  of  the  vehicle,  in  this  case  the  nose  radius,  Rn.  In  fact, 
it  has  been  shown,  e.g. ,  Refs.  1,  2,  20,  25,  that 

- *  O(f)  (i) 

rb 

in  tne  stagnation  region  of  a  hypersonic  vehicle,  ^ased  on  this  assumption,  the  order 

of  magnitude  of  each  term  in  the  expanded  governing  equations  may  be  calculated. 

Retaining  all  terms  from  order  one  to  those  of  order  (Re  )_1  yields  the  equations 

s 

appropriate  for  a  low  Reynolds  number  analysis.  Utilizing  the  method  of  separation 
of  variables,  the  four  non-linear  partial  differential  equations  are  reduced  to  a  tenth 
order  set  of  six  non-linear  ordinary  differential  equations  with  split  boundary  condi¬ 
tions.  The  primary  additional  assumptions  are:  1)  the  shock  is  concentric  with  the 
body,  and  2)  the  variables  are  separable.  The  resulting  system  of  equations  and 
boundary  conditions  is  then  utilized  to  obtain  first  order  non-similar  solutions  to  the 
flow  field  in  the  forward  region  of  the  hypersonic  low  Reynolds  number  viscous  layer. 
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II.  BASIC  RELATIONS 


The  utility  of  the  macroscopic  conservation  equations  which  lead  to  the 
boundary  layer  equations  in  the  high  Reynolds  number  continuum  regime  has  been 
amply  demonstrated  over  the  past  fifty  years.  In  extending  the  analysis  to  lower 
Reynolds  numbers,  it  is  assumed  that  the  density  is  sufficiently  large  that  the  same 
initial  macroscopic  governing  equations  are  applicable.  In  the  absence  of  external 
force  fields  the  steady  state  form  of  these  equations  is: 


Conservation  of  Species  i 

V  •  (pjv.)  =  w4 

(2) 

Conservation  of  Momentum 

•  v)“v  =  -  vp  +  y  •  r 

(3) 

Conservation  of  Energy 


~v  ■  v  . 

V  (h  +  — — )  =  -  V 


Q  +  v  •  (r  •  v) 


(4) 


Equation  of  State 

P  =  P  (h,  P)  (5) 

Summation  over  all  species  in  equation  (2)  yields  the  global  continuity  equation 
V  •  (p'T)  =  0  (6) 

Implicit  in  the  above  definitions  is  that  the  gaseous  system  is  composed  of  a  mixture 
of  chemically  reacting  perfect  gases.  Upon  neglecting  the  Dufour  effect  and  the  trans- 
portof  radiative  energy,  the  heat  flux  vector  Q  is  (Ref.  26), 

Tf=  -  kvT  +  Lp?V.h.  (7) 

i*i  li 


i.  e. ,  the  sum  of  the  Fourier  conduction  and  the  energy  transferred  by  diffusion.  From 
the  definition  of  h,  it  is  seen  that 

vT  =  vh  -  Eh.vC.  (8) 

n  *-  i  J 


where  is  the  frozen  specific  heat: 
P  __ 

C  =  Z  C  C 
P  i  i  Pi 


(9) 


As  a  first  approximation,  it  is  assumed  that  the  gaseous  components  with  similar 
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thermal  and  chemical  behavior  may  be  treated  as  a  single  averaged  component  so  that 
the  sum  total  can  be  treated  as  a  binary  mixture  of  "air  atoms"  and  "air  molecules". 
Then,  neglecting  pressure  and  thermal  diffusion  effects,  the  mass  flux  with  respect 
to  the  mass -averaged  velocity  reduces  to  the  classical  Fick's  diffusion  relation 

Ajvci  (io> 

Substitution  of  eqs.  (8)  and  (10)  into  (7)  yields 

=  -  -jjjr  [?h  +  (Le-l)  ShjVCj]  (11) 

Assuming  a  Lewis  number  of  unity  thus  yields 

Q*=-  l£-  vh  (12) 

and  the  diffusion  equations(2)  are  uncoupled,  even  for  the  case  of  non -equilibrium  flow, 
without  resorting  to  any  specific  constraint  on  the  chemical  source  term,  i.  e. ,  the 
term  L^h.  is  accounted  for  implicitly. 

With  these  simplifying  assumptions  the  governing  equations  valid  in  the  low 
Reynolds  number  regime  appear  expanded  below  for  the  coordinate  system  shown  in 
figure  2. 

Continuity: 


(purj)  +  "  <pvrj)  +  KBrjpv  =  0 

x-Component  of  Momentum: 

[du  du  .,  1  dP  d  T  dul 

u -  +  v -  +  K  uv  =  -  — +  -  <j  - 1 

ax  ay  B  J  ax  dy  ByJ 

+  Kb  [(2+i)  >*  ay"  ("u)] 

y-Component  of  Momentum: 


T  av 

,  av 

-  K  u2l  =  - 

JJ2.  +  i_ 

a 

I>1 

“L  ax 

ay 

B  J 

ay 

3 

ay 

ay  J 

+^_f 

>1 

LL  au 

1  r  ay 

2_  _a_ 

r  ( 

£*<?-)] 

ax  [ 

syJ 

3 

ay 

(13) 


(14) 


(15) 
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Energy: 


It  should  be  noted  that  the  governing  equations  presented  here  are  more  com¬ 
plete  than  those  given,  for  example,  in  Ref.  20,  in  that  all  terms  ranging  from  order 

unity  to  order  (Re  )  1  appear.  Although  Ho  and  Probstein  state  that  they  have  consi- 

S  -1 
dered  all  terms  ranging  from  order  unity  to  order  (cRe  )  ,  actually  one  term  of 

-1  s 
order  (*Re  )  seems  to  be  missing  in  their  energy  equation.  However,  this  intro- 
s 

duces  no  error  in  their  computations  since  the  missing  term  vanishes  identically  at 
the  stagnation  point. 

State: 

P  ~  p(h*  p)  (17) 

Viscosity  Law: 

P)  U8) 

At  this  point  then,  there  are  six  unknowns  u,  v,  h,  p,  p  and  ^  which  may  be 
determined  by  solving  the  six  equations,  (13)  through  (18),  simultaneously. 
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III.  BOUNDARY  CONDITIONS 


The  location  of  the  edge  of  the  viscous  layer  within  the  shock  layer,  and  the 
magnitude  of  the  physical  variables  at  the  outer  edge  of  the  viscous  layer  are  not 
known  "a  priori".  Therefore,  for  convenience,  it  was  decided  to  integrate  the  govern¬ 
ing  equations  from  the  wall  to  the  shock  wave,  where  the  physical  variables  can  be 
calculated. 

It  was  assumed  that  the  shock  in  the  stagnation  regi  n  is  a  surface  of  dis¬ 
continuity  concentric  with  the  body.  It  was  further  assumed  that  thermochemical 
equilibrium  is  achieved  behind  the  shock.  Use  of  the  Rankine-Hugoniot  relations 
yield  the  outer  boundary  conditions  behind  the  shock  wave: 


u(x, 

«s>  -  Vv„ 

(19) 

V(x, 

-rKB2x2)<v» 

(20) 

P(X, 

6s)M1-kb2x2)  <1-0P„V? 

(21) 

„  2 

v 

9  9  9 

h(x. 

Z  Z  Z  oo 

6s)  Ml  -  KB  x  )  d  -  O  2 

(22) 

where 

c  = 

*s 

(23) 

It  was  also  assumed  that  there  is  no  slip  or  temperature  jump  at  the  wall 

and  hence  the  boundary  conditions  at  the  wall  are: 

u(x, 

0)  =  u  -  0 
w 

(24) 

v(x. 

0)  =  v 

w 

(25) 

h(x. 

0)  =  h 

w 

(26) 

Note  that  the  number  of  boundary  conditions,  namely  seven, 

equals  the  order  of  the 

mathematical  system. 

Although  the  magnitude  of  the  physical  variables  can  be  calculated  behind 
the  shock,  the  distance  from  the  wall  to  the  shock  is  not  known.  However,  one  may 
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introduce  a  constraint  on  the  family  of  solutions  based  on  the  conservation  of  mass 
which  takes  the  form: 


O  v  r 
w  w 


1+j 


pu(2r)jdy 


(27) 


This  equation  which  serves  to  delimit  the  thickness  of  the  shock  layer  6g  is  based  on  the  fact 
that  the  gas  flow  between  the  body  and  the  shock  wave  is  the  sum  of  the  flow  which  enters 
the  shock  layer  from  the  upstream  side  of  the  shock  wave  plus  the  gas  injected  at  the 
surface. 
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IV.  SEPARATION  OF  VARIABLES 


The  governing  partial  differential  equations  may  be  solved  by  means  of  the 
technique  known  as  separation  of  the  variables,  which  has  previously  proven  useful 
in  treating  the  low  Reynolds  number  viscous  layer,  (Refs.  1,  17,  20,  21).  In  the 
manner  suggested  by  Probstein  and  Kemp  (Ref.  17)  the  following  separable  forms  of 
the  dependent  variables  were  assumed  which  are  compatible  with  the  boundary  con¬ 


ditions  at  the  shock  wave. 

U  =  KBXUl(y)  <28) 

v  =  (1  -  J-  Kb2x2)  v^y)  (29) 

b  =  (1  -  KB2x2)h1(y)  -  KB2x2h2(y)  (30) 

*>  =  <!  “7  KB2x2)Pl(y)  (31) 

P  =  (1  -  Kb2x2)  Pl(y)  -  KB2x2p2(y)  (32) 

li  -  (1  -  Kb2x2)  Ml(y)  (33) 

r  -  x(l  +  KBy)  (34) 


where  the  terms  with  subscripts  1  or  2  are  functions  of  y  alone.  For  convenience,  the 
following  non-dimensional  forms  are  introduced: 


00 


where  the  subscripts  ®  and  s  refer  respectively  to  the  free  stream  conditions  and  to  the 
state  immediately  behind  the  shock.  Note  that  the  Reynolds  number  is  defined  by: 


KBPs 


(36) 
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For  convenience,  a  compressibility  transformation  is  also  introduced,  i.  c. 

17  ^  Pdy  (37) 

With  the  above  simplification^  the  low  Reynolds  number  viscous  layer  equa¬ 
tions  reduce  to  the  following  system  of  dimensionless  non-linear  ordinary  differential 
equations  with  variable  coefficients  applicable  for  first  order  non-similar  solutions 
in  the  forward  region  of  a  hypersonic  vehicle: 

(pvV  +  (1  +  j)  (-7-  +  v)  =  0  (38) 


*P,  +  P 


u(-y  +  v)  4  pvu'  -2-—  2  =  (J^tT)'  +  (2  +  j)  -  (fu )’ 
cpvv>  +  p'  =  y  e(pp7')’  +  (l+D^u*  •  ”§”  <5^*1 
Po  +  11  <  “T~  +”)  =  0 


(39) 

(40) 

(41) 

(42) 

(43) 


Prpvh;*  IfijSkp  +  (1  +  j)  iff; 

ps^  4  2  [j^^acr')'  4  (i  +  j)  pair  -  (pa2)^] 

*  '  +  (1  +  j)  ~P*  **2  +  27~  [“2  +  €5vu*  - (Kj  4"E2)J 

P=/5,(h1,  Pj)  (44) 

1 3  =  7T(h1»  Pj)  (45) 

where  the  primes  denote  derivatives  with  respect  to  77.  It  is  seen  that  separation  of  the 
variables  has  increased  the  number  of  unknowns  to  eight,  including  u,  V,  TT  ,  T  , 
p  ,  P9.  p»  p.  This  requires  that  the  number  of  governing  equations  also  be  eight, 
equations  (38)  to  (45),  while  the  overall  order  of  the  mathematical  system  is  ten 
requiring  ten  boundary  conditions. 

The  boundary  conditions  become: 

At  the  surface,  77  =  0,  equations  (24),  (25)  and  (26)  yield 

u  =  0,  v  =  v  ,  (46) 

h^  =  h^  =  -  h2  (isoenthalpic  wall) 
w  w 
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At  77  =  77  .  equations  (19)  through  (22)  yield: 
s 

»  -1.  V  --1,7  -O.-Pj  »(l-<),?2  =0  (47) 

S  S  8  S 

The  constraint  equation  becomes: 

a  +  V1  + 1  +  pwvw  -  f- [2d +  y>]  dv  <«) 

Here  it  is  noted  that  for  the  special  case  of  zero  mass  transfer,  i.  e. ,  a  v  =0, 

w  * 

examination  of  eq.  (38)  indicates  that  the  constraint 

“v»  =  0  (49) 

w 

becomes  identical  with  eq.  (48).  This  latter  result  (Ref.  1)  follows  from  the  fact 

that  since  *u  and  V  are  both  zero  (for  zero  mass  transfer),  unlessT1  also  vanishes, 
w  w  w 

one  obtains  the  physically  untenable  result  of  an  infinite  density  gradient  at  the  wall, 
(Ref.  20). 

In  general,  the  transport  properties  depend  on  both  temperature  and  com¬ 
position.  However,  since  the  diffusion  equations  have  been  uncoupled,  the  compo¬ 
sition  of  the  gas  is  not  calculated  explicitly,  and  so  it  is  necessary  to  resort  to  further 
approximations.  Constant  Prandtl  and  Lewis  numbers  have  already  been  introduced 
and  hence  the  thermal  conductivity  and  the  diffusion  coefficients  are  not  required 
explicitly.  Thus,  only  the  viscosity  coefficient  remains  to  be  evaluated. 

It  is  assumed  that  the  gas  is  in  chemical  equi’ibrium.  or  more  explicitly 
that  the  state  of  the  gas  can  be  determined  from  a  Mollier  diagram  (e.g. ,  Ref.  27). 
That  is,  knowledge  of  two  state  functions  (e.g. ,  pressure  and  enthalpy)  uniquely  deter¬ 
mines  the  others  (e.g.,  temperature,  composition,  density,  etc.).  As  already  noted 
it  is  assumed  here  that  the  viscosity  coefficient  can  be  approximated  by  Sutherland's 
formula  for  air  (Ref.  28) 

M  '  Md)  -  1. 16  x  l(f5  (-JV)  (£-)  3/2  (50) 

with  T  in  R,  ^  is  given  in  lb.  /ft.  sec.  However,  since  (from  the  Mollier  diagram) 

T  =  T(h,  p)  (51) 
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then  implicitly 


p  =  P)  (52) 

and  the  equation  of  state  may  also  be  obtained  from  the  Mollier  diagram,  i.e. , 

P  =  p(h,  p)  (53) 

Conditions  behind  the  shock  were  determined  from  shock  tables  which  in¬ 
corporate  the  1959  ARDC  model  atmosphere,  (Ref.  29). 
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V.  DISCUSSION  OF  RESULTS 


Solutions  were  obtained  on  an  IBM  7094  high-speed  digital  computer  in  the 

hypersonic  low  Reynolds  number  viscous  layer  flight  regime.  Calculations  were 

carried  out  for  a  range  of  flight  speeds  from  15,  000  ft.  /sec.  to  25,  000  ft.  /sec. ,  an 

altitude  range  from  150,  000  ft.  to  350, 000  ft.  for  Reynolds  numbers  which  were 

varied  between  25  and  15,000  and  a  mass  transfer  parameter,  -0.4  £  f  s  0.  These 

w 

results  were  then  compared  with  earlier  theoretical  studies  in  the  continuum  regime 
and  in  the  low  Reynolds  number  viscous  layer  regime. 

The  non-dimensional  system  of  equations,  (38)  to  (45),  was  integrated  from 
the  wall  out  to  the  shock  wave.  Since  the  boundary  conditions  are  split,  it  is  necessary 
to  employ  an  iterative  procedure.  Some  typical  solutions  for  the  viscous  shock  layer 
profiles  of  “u,  "v,  IT  and  p”  are  shown  in  Figures  3  to  6.  As  indicated  by  previous  stu¬ 
dies,  the  viscous  effects  extend  further  into  the  shock  layer  with  decreasing  Reynolds 
number. 

After  numerous  calculations  for  a  wide  range  of  altitudes,  flight  speeds, 
Reynolds  numbers,  wall  temperature,  and  mass  transfer  rates,  it  was  found  that 
all  of  the  data  could  be  correlated  with  several  standard  non-dimensional  ratios  as 
functions  of  the  Reynolds  number  in  combination  with  the  shock  density  ratio. 

All  of  the  results  of  this  study  were  normalized  with  respect  to  boundary  layer 
predictions  for  precisely  the  same  flight  conditions  and  gas  model.  New  boundary 
layer  solutions  for  this  study  were  obtained  by  retaining  only  first  order  terms  in 
equations  (38)  to  (45),  thereby  yielding  the  high  Reynolds  number  asymptote.  These 
boundary  layer  results  are  in  reasonable  agreement  with  those  of  other  investigators, 
(Ref.  33).  However,  the  discrepancies  between  the  various  predictions  is  nearly  as 
large  as  the  low  Reynolds  number  effects.  Thus,  in  order  to  assess  the  effects  of  low 
Reynolds  number  alone, we  obtained  these  new  boundary  layer  solutions  which  are  com¬ 
pletely  compatible  with  our  low  Reynolds  number  solutions.  The  details  of  these  calcu¬ 
lations  will  appear  in  a  paper  to  be  presented  shortly  in  which  the  hypersonic  shock 
layer  is  examined  over  a  region  in  which  continuum  concepts  may  reasonably  be  ex¬ 
pected  to  be  valid,  i.  e. ,  from  boundary  layer  concepts  to  merging  of  the  shock  with 
the  fully  viscous  flow  behind  it. 
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All  of  the  normalized  heat  transfer  rates  with  no  mass  transfer  obtained  for 

the  present  investigation  are  shown  in  figure  7.  Note  that  two  values  of  the  heat 

2 

transfer  rate  are  given  for  a  particular  value  of  e  Reg.  This  is  due  to  the  assumed 
form  of  the  mathematical  separation  of  variables.  Combining  the  two  algebraically 
yields  a  first-order  non-similar  solution  to  the  heat  transfer  rate  partly  around  the 


and  is  valid  for 

sin2(Kfix)  s  0(c) 


(55) 


For  example,  at  an  altitude  of  150, 000  ft. ,  a  flight  speed  of  15,  000  ft.  /sec. ,  c  %  0(0. 1) 
and  so  (Kgx)max  18. 5°.  It  is  seen  that  the  concept  of  similarity  is  a  reasonably  good 
approximation  at  higher  Reynolds  numbers,  in  agreement  with  Ref.  34.  However,  when 
the  viscous  effect  are  significant  over  a  major  portion  of  the  shock  layer,  this  assump¬ 
tion  becomes  less  reasonable,  which  is  not  in  agreement  with  Ref.  34.  However,  since 
the  results  of  the  present  study  are  only  valid  in  the  larger  Reynolds  number  range  of 
the  low  Reynolds  number  regime  (i.  e. ,  viscous  layer  subregime)  not  much  else  can  be 
said  here  concerning  the  discrepancy. 

It  was  found  that  plotting  all  of  our  normalized  heat  transfer  rate  results 

2 

with  respect  to  the  parameter  c  Reg  reduced  the  scatter  to  a  minimum.  Furthermore, 
it  reduced  the  discrepancies  between  various  theoretical  and  experimental  studies  pre¬ 
viously  conducted  (Refs.  1,  20,  23,  30,  31,  35  and  36).  It  is  seen  that  the  results  of  the 
present  analysis  indicate  the  same  general  trends  but  a  somewhat  larger  magnitude  than 
all  of  these  earlier  studies  with  the  exception  of  that  of  Ferri,  et.  al.  (Ref.  30)  which 
is  almost  identical  with  our  present  results.  Coincidently,  only  in  Ref.  30  and  in  the 
present  study  was  a  self  consistant  boundary  layer  asymptote  determined  in  order  to 
separate  the  low  Reynolds  number  effects.  However,  it  should  also  be  noted  that  the 
results  of  Ref,  30  do  not  increase  with  decreasing  Reynolds  number  as  rapidly  as  the 
present  more  complete  analysis  indicating  the  need  to  retain  more  higher  order  terms 
with  decreasing  Reynolds  numbers.  This  is  further  demonstrated  by  the  results  of 
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Kao  (Ref.  35),  i.e. ,  the  divergence  of  his  second-order  boundary  layer  and  third- 

order  boundary  layer  predictions  at  small  Reynolds  numbers,  and  comparison  with 

Ho  and  Probstein’s  results  (Ref.  20).  It  is  noted  that  the  conditions,  gas  model  and 

boundary  layer  asymptote  arc  the  same  in  both  of  these  references. 

It  is  stressed  that  the  calculations  in  this  study  were  carried  out  down  to 
2  3 

values  of  Re  ,  c  Re  and  c  Re  below  the  range  in  which  it  is  expected  that  the 
s  s  s 

viscous  layer  will  merge  with  the  shock  wave,  in  order  to  distinguish  between  the 

various  predictions.  From  Refs.  20  and  30,  it  would  appear  that  the  lower  limit 

of  the  viscous  layer  regims  is  of  the  order  of  Re  =  100.  Ho  and  Probstein  (Ref.  10) 

8 

extended  their  viscous  layer  results  into  the  merged  viscous  layer  regime  by  modi¬ 
fying  a  previous  study  (Ref.  17,  constant  density  solution).  In  a  more  recent  study 
by  Levinsky  and  Yoshihara  (Ref.  21),  the  same  system  of  equations  and  perfect  gas  model 
employed  by  Probstein  but  utilizing  a  different  gas,  was  solved  from  the  free  stream  through 
the  merged  shock  and  viscous  layer  to  the  wall.  Their  results  are  naturally  similar,  except 
that  paradoxically  they  are  shifted  by  more  than  an  order  of  magnitude  to  higher  Reynolds 
numbers  (see  Ref.  1).  However,  the  peak  in  their  results  with  a  downward  trend  at 
the  lower  Reynolds  numbers  indicates  the  merged  shock  layer  structure.  As  noted 
earlier,  Cheng  (Ref.  22)  did  not  make  the  poor  assumption  of  a  perfect  gas  shock 
structure,  but  modified  the  Rankine-Hugoniot  shock  relations  to  account  for  the 
normal  components  of  the  energy  and  momentum  transfer  through  the  shock  wave. 

His  results  are  also  shown  in  figure  7. 

It  is,  of  course,  expected  that  a  realistic  merged  viscous  layer  analysis 

will  yield  results  which  will  show  a  downward  trend,  somewhat  similar  to  Cheng's, 

2 

and  will  not  continue  to  increase  with  decrease  in  *  Re  as  the  various  viscous 

s 

layer  models  indicate  in  figure  7. 

Recent  experimental  results  obtained  by  Hacker  and  Wilson  (Ref.  31)  are  also 
shown  in  figure  7  utilizing  their  correlation  relation 


Q 


(Q), 


=  l  + 


0.4  ±  0.1 


(56) 


rBL  V  Res 

and  assuming  e  =  0. 12.  Only  minor  low  Reynolds  number  effects  are  indicated, 
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whereas  Ferri,  et.  ai.  *s  (Ref.  30)  experimental  results  appear  to  justify  their  own 
theory  indicating  a  more  significant  low  Reynolds  number  effect.  The  more  recent 
experimental  results  of  Vidal  and  Wittliff  (Ref.  36)  also  agree  well  with  those  in 
Ref.  30  iidicating  the  more  significant  low  Reynolds  number  effect. 

The  numerical  discrepancies  between  the  results  of  the  present  study  and  those 
obtained  earlier  by  us  (Ref.  1)  is  mainly  due  to  the  method  of  calculation.  The  results 
in  Ref.  1  were  obtained  on  an  analog  computer,  whereas,  the  present  results  were 
obtained  in  a  more  exact  fashion  on  an  IBM  7094  digital  computer  which  is  inherently 
more  accurate  than  any  analog  computer. 

Normalized  skin  friction  results  of  this  investigation  are  compared  with  some 
others  reported  on  in  the  open  literature  in  figure  8.  It  is  seen  in  this  figure  that  the 
present  study  shows  a  larger  low  Reynolds  number  effect  than  earlier  studies. 

It  was  shown  in  Ref.  1  that  the  gas  model  employed  was  naturally  not  as  signi¬ 
ficant  in  the  calculations  of  skin  friction  as  for  heat  transfer  rates.  In  fact,  the  low 
Reynolds  number  nature  of  the  flow  was  shown  to  be  the  important  criteria.  Most 
of  the  terms  neglected  in  the  studies  of  Ho  and  Probstein  (Ref.  20)  appear  in  the  equa¬ 
tion  for  the  x-component  of  momentum,  which  is  most  influential  in  the  determination 
of  the  shear  stress.  This  can  also  be  seen  by  the  discrepancies  in  the  two  predictions 

3 

of  Kao  (Ref.  35)  and  Ho  and  Probstein  (Ref.  20).  It  found  that  the  abscissa  e  Re 

s 

used  in  this  figure  reduced  the  scatter  such  that  most  of  the  calculations  were  right  on 
the  curve  and  the  rest  had  a  scatter  of  no  more  than  two -three  percent. 

It  has  already  been  shown  that  in  the  absence  of  mass  transfer  it  is  convenient 
to  normalize  the  heat  transfer  rate  and  surface  shear  stress  by  means  of  the  high 
Reynolds  number  (boundary  layer)  asymptote.  This  suggests  that  when  mass  transfer 
is  present  one  can  again  normalize  the  low  Reynolds  number  results  by  means  of  the 
corresponding  high  Reynolds  number  asymptote.  A  summary  of  the  effects  of  air  injection 
into  air  at  high  Reynolds  numbers  from  new  boundary  layer  calculations  is  shown  in 
figure  9,  where  the  mass  transfer  parameter  B  is  defined  as 


B  = 


(H  -h  ) 
s  w' 


(57) 


18 


In  figures  10  and  11  the  normalized  (i.  e. ,  normalized  with  respect  to  our  new 

boundary  layer  solutions)  contributions  of  the  heat  transfer  are  shown  as  a  function  of 
2 

c 'Re  ,  wall  temperature  and  dimensionless  mass  transfer  rate  f  defined  a  ; 
s  w 


f  =  - 
w 


0  v 

^w  w 


V  (l  + « *\Av  ~±r 

It  is  seen  that  the  effectiveness  of  the  mass  transfer  process  continues  into  the  low 
Reynolds  number  regime,  being  somewhat  less  effective  for  a  cool  wall. 

Interestingly  the  single  curve  for  the  present  results  of  normalized  skin 
friction  shown  in  figure  9  represents  not  only  the  zero  mass  transfer  cases  but  also 
all  of  our  results  with  mass  transfer. 


As  mentioned  earlier,  calculations  at  Reynolds  numbers  as  high  as  15,000 
were  carried  out  for  this  study.  Our  first  attempts  at  normalization  (using  the  boundary 
layer  correlations  of  Scala  and  Gilbert,  Ref.  37)  indicated  that  at  this  high  a  Reynolds 
number  the  low  Reynolds  numbers  effects  on  heat  transfer  and  skin  friction  were  quite 
small.  However,  there  were  discrepancies  such  that  satisfactory  correlations  could 
not  be  achieved.  Thus,  it  was  decided  to  carry  out  a  compatible  boundary  layer 
analysis.  Using  these  new  solutions  we  have  been  able  to  correlate  all  of  the  low 
Reynolds  number  results.  However,  it  is  now  seen  that  even  at  this  high  a  Reynolds 
number  that  there  is  a  difference  in  the  zero  mass  transfer  heat  transfer  rate  pre¬ 
dictions  of  about  ten  percent  and  about  thirty  percent  for  the  skin  friction.  It  is,  thus, 
seen  that  the  low  Reynolds  number  effects  appear  to  extend  to  much  higher  Reynolds  numbers 
than  had  previously  been  reported.  This  will  be  discussed  in  greater  detail  in  another 
paper  soon  to  be  released. 

In  the  high  Reynolds  number  regime,  the  shock  detachment  distance  is  a  function 
of  body  shape,  density  and  velocity,  independent  of  surface  temperature  and  mass 
transfer  at  moderate  rates  of  injection.  However,  at  lower  Reynolds  numbers,  the 
shock  detachment  distance  is  a  function  of  body  shape,  flight  speed,  altitude,  Reynolds 
number,  wall  temperature  and  mass  transfer  rate.  This  effect  was  investigated  in  order 
to  predict  the  effects  of  mass  transfer  in  the  low  Reynolds  number  viscous  layer  regime 
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upon  shock  wave  detachment  distance.  Figure  12  shows  the  variation  of  6  /R  with 

S  13 

the  above  parameters  at  a  single  flight  speed  and  two  wall  temperatures.  Intuitively, 

one  expects  that  due  to  the  extended  interaction  region,  mass  injection  will  push  the 

shock  wave  further  away  from  the  body  with  decreasing  Reynolds  numbers.  As  shown 

in  figure  12,  this  is  indeed  found  to  be  the  case. 

As  noted  earlier,  the  location  of  the  "edge”  of  the  boundary  layer  within  the 

shock  layer  is  not  known  Ma  priori".  However,  "a  posteriori"  the  limit  of  the  viscous 

effects  6«t  /R^  can  be  defined  as  the  location  at  which  o  /o  =  0. 995.  Using  this 

definition,  6  /R  for  the  same  conditions  is  also  shown  in  figure  12.  The  terminus 
BL  B 

at  the  lowest  Reynolds  number  is  the  intersection  of  the  edge  of  the  boundary  layer 

with  the  shock  detachment  distance.  Obviously,  this  theory  fails  at  a  somewhat  higher 

Reynolds  number  since  the  shock  wave  has  also  thickened  appreciably  before  merging 

with  the  viscous  layer.  Also  shown  is  the  estimate  based  on  our  new  boundary  layer 

solutions.  It  is  seen  that  even  at  the  large  value  Re  =  15,  000  the  boundary  layer  thick- 

s 

ness  prediction  is  greater  for  boundary  layer  analysis. 
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VI.  CONCLUSIONS 


1)  Useful  correlations  for  heat  transfer  and  skin  friction  have  been  obtained 
for  the  hypersonic  low  Reynolds  number  viscous  layer. 

2)  The  effects  of  mass  transfer  of  air  molecules  into  the  viscous  layer 
have  also  been  determined  and  correlated  herein. 

3)  The  range  of  applicability  of  the  present  model  has  been  delimited. 
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Figure  1.  Hypersonic  Flight  Regimes 


Figure  2.  Coordinate  System 


Figure  3.  Profiles  Re  =  15,000 
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Figure  4.  Profiles  Re 


Figure  6.  Profiles  Re 


1.  PRESENT  CORRELATION  FOR  0 

2.  PRESENT  CORRELATION  FOR  0 

3.  ANALOG  SOLUTION  (REF.  0 


Figure  7.  Normalized  Heat  Transfer 


Figure  8.  Normalized  Skin  Friction 


Figure  9.  Boundary  Layer  Mass  Transfer  (Q  /Q 


Figure  10.  (Q  )/(Q  )__  vs.  €  Re  (Mass  Transfer) 
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Figure  1Z.  Shock  Detachment  and  Viscous  Layer  vs  Re 
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